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Abstract 

It is well known that in low energy supersymmetry the ratio of the gaugino mass to the 
gauge coupling squared, M/a, is renormalization group invariant to one-loop. We present a 
systematic analysis of the corrections to this ratio, including standard two-loop corrections 
from gauge and Yukawa couplings, corrections due to an additional U{1)' gaugino, threshold 
corrections, superoblique corrections, corrections due to extra matter, GUT and Planck scale 
corrections, and "corrections" from messenger sectors with supersymmetry breaking commu- 
nicated via gauge-mediation. We show that many of these effects induce corrections at the 
level of a few to tens of percent, but some could give much larger corrections, drastically 
disrupting the renormalization group extrapolation of the ratio to higher scales. Our analysis 
is essentially model- independent, and therefore can be used to determine the ambiguities in 
extrapolating the ratio in any given model between the weak scale and higher scales. 



1 Introduction 



It is well known that any hope of stabilizing a hierarchy of scales requires super symmetry |jl|. 
The disparity between the weak scale and the grand unified theory (GUT) scale or the Planck 
scale (hereafter lumped together as the "high scale" ) is bridged in a softly broken super symmetric 
standard model by rendering scalar mass renormalization to be at most logarithmically diver- 
gent. The canonical approach to connecting these disparate scales is through the renormalization 
group (RG), and there has been enormous effort in calculating and evaluating the RG equations 
analytically and numerically, e.g. Refs. |, |, |, |, |, |, |, |, 0, 0, Most of this effort 
has been directed towards either the consequences of a given high scale theory (be it a GUT, 
string theory, or other high scale proposal) on weak scale phenomenology, or the consequences 
of present (or proposed) weak scale measurements on high scale physics. 

A general, softly broken, weak scale supersymmetric model has a large number of additional 
parameters beyond those of the standard model (SM) Most of the parameters are inter- 
connected through the RG equations, and therefore successful extrapolation generally requires a 
simplifying framework (or an organizing principle) to be imposed either at the weak scale or the 
high scale. Even with the varied simplified frameworks that have been traditionally used, such as 
supergravity-motivated models, gauge-mediated models, and so on, there are a sufficient number 
of parameters that one usually resorts to numerical sampling. While this approach usually gives 
the correct general trends, it could also easily give misleading results if the physics of the model 
depends very sensitively on a sampled parameter. 

Disentangling the dependencies of any given supersymmetric parameter on other parameters 
has proved to be formidable task. The first step, extracting weak scale observables, is likely going 
to be challenging both experimentally (discovery, and then mass and coupling measurements) 
and theoretically (mixings between gauginos, one-loop corrections, etc.). Hence, the subsequent 
extrapolation of the full theory to higher scales will be unreliable unless quantities can be found 
that are not highly interdependent on other parameters of the theory. One such quantity is the 
ratio of the gaugino mass to the gauge coupling squared, Ma/aa, where a = 1, 2, 3 corresponds to 
J7(l)y, SU{2)l, SU{3)c- In principle, one can extract Ma/aa by "merely" measuring the gaugino 
mass, given that the standard model gauge couplings are by now well measured. Remarkably, 
the one-loop RG equation for Ma/oa vanishes |13], which suggests that a measurement at the 
weak scale would determine the value at the high scale, up to small two-loop corrections (albeit 
scaled by ~ In Qhigh/Qweak)- 

The usefulness of the ratio Ma/oa at higher scales is, ultimately, model-dependent. However, 
there are strong motivations to think that the pattern of the gauge couplings and the gaugino 
masses at the high scale will determine a great deal about the high scale theory. In GUTs, 
one ordinarily expects that the gauge couplings and the gaugino masses will "unify", or take 
the same value, after being embedded in the GUT group [14|. Even in string theory there are 
reasons to suggest the gaugino masses could be unified in the simplest scenarios If we could 
reliably calculate Ma/aa at the high scale, no doubt one could match string phenomenology to 
such quantities, strengthening the predictions of successful models (and eliminating classes of 
unsuccessful models). In addition, given a motivation for thinking i^gut or (^string should take 



2 



a certain value, it would also be trivial to extract the gaugino masses by themselves from the 
ratios Ma/aa- 

The difficulty of unambiguously extrapolating the individual quantities ga and M^, or even 
other weak scale parameters such as squark and slepton soft masses, is mainly due to the one- 
loop RG evolution and one-loop dependence on threshold corrections. In particular, extrapolating 
up 13 orders of magnitude in scale to the apparent unification scale Munif ~ 2 x 10^^ GeV is 
notoriously complicated by weak scale supersymmetry threshold corrections, GUT scale threshold 
corrections, extra matter at intermediate scales, Planck scale corrections, etc. Using the ratio 
Ma/aa mitigates many of these issues due to the one-loop RG invariance, as we will see. 

Our goal is ultimately to determine how well Ma/aa can be known at higher scales, and 
thus we will be primarily concerned with running up to the high scale (or messenger scale), 
rather than the more traditional approach of running down from a boundary condition imposed 
by some supersymmetry breaking ansatz. Our results can, of course, be easily inverted should 
there be a motivation to do so. Bottom-up approaches to supersymmetry have been considered 
before (see e.g. 0, 16, 17, |l^]), generally with the philosophy that weak scale phenomenology 



ought to be extrapolated to the high scale, leaving the high scale matching to GUTs or string 
theory. Our approach is similar, in that we start from the weak scale, but differs from some of 
the previous analyses by extrapolating only the ratios Ma/aa, and yet also considering various 
high scale effects. Throughout most of this work we illustrate the results of extrapolation from 
near the weak scale ~ 1 TeV up to the apparent unification scale ~ 10^^ GeV; should the soft 
mass generation scale be much lower as in gauge-mediated supersymmetry breaking models, the 
size of the corrections will be smaller, but can still be easily extracted from most of the graphs we 
present. Similarly, effects in high scale theories (with gravity-mediated supersymmetry breaking) 
can give corrections to the ordinary expectations for Ma/aa at the weak scale. By quantifying 
the various effects that could disrupt the one-loop expectations for Ma/aa, we are simultaneously 
showing that a naive interpretation of the gaugino mass ratios at the weak scale could easily give 
misleading conclusions, but a systematic evaluation of the effects discussed here could lead to 
intriguing signals of physics at much higher scales. 

The most obvious class of corrections to the one-loop invariant Ma/aa are the two-loop 



corrections |19, 20, 21 1 (in the super symmetric DR renormalization scheme ^^). The two- loop 
corrections ought to be small, if perturbation theory is validQ. However, there are reasons to 
suggest that the size of two-loop corrections could be larger than one might naively guess. The 
two-loop coefficients are not of 0(1), but in fact 0(10) (in the minimal supersymmetric standard 
model (MSSM) the two-loop coefficients lie in the range 9/5 < i?^^ < 25). Also, two-loop terms 
proportional to the Yukawa couplings are present, and of course the top Yukawa is large at the 
weak scale. Finally, extrapolation to the high scale involves running over 13 orders of magnitude 
in scale, and therefore the (resummed) logarithm from RG evolution is comparatively large. 
(Analogously, if the soft mass generation scale is much lower, then the effect of RG evolution is 
not nearly as large.) 

If two- loop corrections become important, then three- loop corrections may also be important, 
if for no other reason than to check that perturbation theory is indeed valid. However, a complete 

^ Which we assume throughout this paper. 
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Table 1: Dynamical parameter dependence in supersymmetry for one-, two-, and three-loop 
renormalization group equations. In the MSSM the labels i,a = 1,2,3 specify gauge couplings 
and w,x = u,d,e specify Yukawa couplings, with i and w fixed while a and x are summed over. 



three-loop analysis would entail using three- loop gauge and gaugino mass /5-functions, at least 
two-loop Yukawa couplings and scalar trilinear couplings, and two-loop thresholds. The three- 



loop gauge coupling /3-functions for a general MSSM were presented in Ref. |23|, and the three- 



loop gaugino /^-functions were recently presented in Ref. u4. Two- loop /3-functions for the 



Yukawa couplings and scalar trilinear couplings are by now well known [£5, 26]. To the best of our 
knowledge two-loop threshold corrections (in a nondecoupling scheme) have not been calculated, 
although one could implement a decoupling scheme whereby the two-loop coefficients are changed 
as each threshold is crossed; this procedure is an approximation that should be straightforward 
to calculate from generalized two-loop /3-functions, although we decline to present it here. An 
additional well-known complexity in going to three-loops is that the /^-functions are no longer 
scheme independent |^^. In any case, we will not attempt to calculate corrections beyond two- 
loop, instead relying on others' calculations ^] of three- loop corrections to guide us in those 
cases where we suspect perturbation theory may be in trouble. 

The central reason why three-loop corrections should be unnecessary (unless a coupling gets 
large) is that there is no dependence on any additional parameters beyond those needed for the 
two-loop corrections. The intriguing pattern of the supersymmetric parameter interdependence 
is given in Table ||, where we show the dependencies arising at one-, two-, and three- loop order 
in the RG equations. The results to three-loop have been explicitly calculated for the gauge 
couplings g, Yukawa couplings Y, and gaugino masses M ||2^, and we suspect the three-loop 
result for the scalar trilinear couplings A holds based on two- loop results ||2^, p6[ . 

The parameter dependencies in Table || are central to our analysis. They imply, specifically, 
that the set of parameters {g, M, Y, A) are sufficient to compute RG evolution to three-loop (and 
probably to all orders), and that we need not be concerned with the squark or slepton soft masses 
or the bilinear couplings (such as fi or B) directly^. This interdependence does, however, deserve 
a few further comments. 

In the exact supersymmetric limit, only gauge couplings and superpotential parameters 
(including Yukawa couplings) remain. It is therefore not surprising that the gauge couplings and 
Yukawa couplings depend on each other, but not on any soft breaking parameter. That the soft 
breaking parameters M and A depend on only the four parameters {g, M, Y, A) can ultimately 

^The extraction of the scalar trilinear couplings at the weak scale, however, can presumably be obtained only 
by measuring the off-diagonal elements of the sfermion mass matrices. 
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be related to theoretical properties of the superpotential. There has been a significant renewed 
interest in exact /3-functions in (N = 1) super symmetric theories. Using a fully Wilsonian 
treatment, exact one- loop /3-functions for the gauge couplings |29| and gaugino masses 24 1 can 
be found by expressing the soft breaking masses as spurious and then exploiting the holomorphy 
of the Lagrangian. Although these results are not directly useful for calculations with canonically 
normalized fields, they do allow the computation of higher order /3-functions including some of 
the three-loop results mentioned above. 

The outline of the paper is as follows. In Sec. ^, we present two-loop effects in the MSSM. 
We discuss the two- loop RG equations including both the "pure gaugino" and Yukawa terms, 
and the dependencies on the sign of the scalar trilinear coupling as well as the gaugino mass. We 
also discuss weak scale threshold corrections and superoblique corrections to the ratio Ma/aa- 
In Sec. 1^ we discuss a series of effects beyond the MSSM, including an extra Consistently 
adding an extra U{1)' to the gauge structure of the MSSM requires the nontrivial cancellation of 
[/(!)' anomalies, and we relegate of discussion of this to Appendix We consider two models 
(with anomaly cancellation); a minimal U{\)' model and an model. We also consider adding 
extra matter to the theory, and its effect on the ratio Ma/oia- In Sec. Q we discuss high scale 
corrections to the ratio Ma/ da from GUT threshold corrections, effects due to breaking GUTs 
at other scales, Planck scale corrections, and supergravity effects. In Sec. ^ we discuss the 
consequences of extra matter with supersymmetry breaking masses that induce "corrections" via 
gauge-mediation. Finally, in Sec. |6| we present our conclusions. In Appendix we write the 
full two-loop RG equations for the gauge couplings and gaugino masses including all Yukawa 
couplings appropriate to the model, and then apply the results to the MSSM, MSSM + U{1)' 
model, and an model. 



2 Effects in the MSSM 



The first class of corrections to consider are those arising from the "standard" two-loop terms 
in the RG equation for Ma/aa- This will be done by starting with certain restrictive assumptions 
about the initial conditions and the two-loop terms to gain insight into the corrections, and then 
removing the assumptions one by one. Both weak scale threshold corrections and superoblique 
corrections are considered here, since they are generic in any weak scale super symmetric model. 

2.1 Two-loop effects in the MSSM 

To begin and start with a point of reference, consider the RG evolution of Ma/aa using 
the two-loop RG equations in the MSSM. The RG equation for the ratio Ma/aa is given in 
Eq. (^), and throughout much of the following discussion we use "pure gaugino" and "Yukawa" 
contributions to refer to the first and second set of terms in the square brackets of the RG 
equation. First, take the unrealistic case where the Yukawa couplings are set to zero, which will 
allow us to analyze the pure gaugino contributions to the RG equations with the least number 
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Figure 1: Evolution of Ma/aa as a function of scale. The solid, dotted, and dashed lines cor- 
respond to a = 1,2,3 respectively. Note that we have normalized the y-axis by dividing by the 
value Mi/ai at 1 TeV. 



of complexities. Furthermore, let's start by assuming that 

Mi/ai = M2/a2 = Ms/as 



(1) 



holds at a scale Q near the weak scale, but above all weak scale thresholds (these will be discussed 



in Sec. 2^). The evolution of the gauge couplings and gaugino masses is governed by six coupled 
RG equations (see Appendix^), but this is still analytically rather complicated. In Fig. we 
show the evolution of the ratio Mq/qq from 10^ — > 10^^ GeV. By choosing to start the running at 
1 TeV with Eq. (||) as the weak scale boundary condition, we are restricting Mi(l TeV) < 190 
GeV, to prevent the gluino mass Ms^l TeV) from exceeding 1 TeV. For smaller values of Mi 
(and correspondingly M2 and M^), the results are identical since there is only one scale in the 
problem (any one of the gaugino masses determines the scale). Nevertheless, for definiteness we 
take Mi(l TeV) = 100 GeV, M2(l TeV) ^ 190 GeV, and M^il TeV) ^ 527 GeV in the following 
unless otherwise specified. The latter values (M2 and M3) are considered mildly approximate 
since their determination requires the gauge couplings at 1 TeV, and therefore subject to weak 
scale threshold corrections. 

There are a few important remarks to be made based on this simple exercise. First, the size 
of the pure gaugino corrections to the one-loop invariant Ma/ota are less than 5% after scaling 
13 orders of magnitude. In all cases the correction is a slight increase with increasing scale 
(RG equation is positive), reflecting the usual choice of phases in the initial condition, Eq. (|^). 
The most striking feature is that the size of the two-loop pure gaugino corrections are, in order 
from largest to smallest (at any scale), M2/a2) -^i/c^i) ^3/0^3- That such an ordering should 
be expected can be seen by examining the RG equation for the difference of two ratios (again, 
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Yukawa couplings set to zero), 

- 2 / (2) 




{b^^J - (2) 



where the differences are 



d(2) _ n{2) _ /_154 98 32 N 

-°2c -Die — V 25 ' 5 ' 5 ^ 

r(2) _ r{2) _ /144 _18 18 N 

-°lc ^3c — V 25 ' 5 ' 5 / 

using the two-loop coefficients given in Appendix Since the difference between the two-loop 
coefficients is positive for the dominant g^M^ term, the RG equations for the difference of the 
ratios (M2/a2 — Mi/ai) and (Afi/ai — M^/a^) are also positive, and thus M2/a2 obtains the 
largest two-loop correction, followed by Mi/ai and M^/a^. This ordering holds if ^2-^2 ^ 53-^3; 
and it would only be if 52-^2 ^ 53-^-^3 and/or g^Mi » g2M2, g'^M^ that one would expect the 
ordering of the size of the two-loop corrections to be different. 

The next stage is to restore the two-loop terms proportional to the Yukawa couplings. In the 
RG equation for the ratio Ma/aa, each Yukawa coupling (squared) is multiplied by the associated 
soft breaking term, the scalar trilinear coupling, which provides the mass scale. Unlike the "pure 
gaugino" terms, additional scales independent of the gaugino masses can partially determine 
the RG evolution. Potentially, one new scale is introduced for every nonzero Yukawa coupling. 
However, we will assume that the 3x3 Yukawa coupling matrices Y„, 'Y^, and Ye (in flavor space) 
can be reduced to the dominant third generation couplings Yt, Yf,, and Yr, and consequently 
the only relevant scalar trilinear couplings are At, Af,, and Ar- The RG equations for the 



Yukawa couplings and scalar trilinear couplings are given in Appendix A. 3 . The size of the terms 
proportional to the Yukawa couplings depend on the competition between terms proportional 
g'^M versus those proportional to Y^A. It is therefore important to recognize that the two-loop 
corrections to Ma/ua need not be proportional to (or suppressed by) a gauge coupling squared, 
but instead a Yukawa coupling squared. 

The Yukawa couplings are extracted at the weak scale from the fermion mass m{m), but they 
also depend on tan /3, the ratio of vacuum expectation values of the two neutral Higgs doublet 
scalar fields. The perturbative lower and upper limits on tan /3 can be found by successfully 
extrapolating the Yukawa couplings up to the high scale Qhigh without encountering a Landau 
pole. We will be considering three cases in the following: one case each of low and high tan (3, 
and in subsequent discussion we usually take an intermediate tan /? value. 

In the low and intermediate tan (3 cases, we can safely neglect the terms proportional to Yi, 
and Yr- If we once again use the initial condition Eq. (0), normalizing to M^/a^ at Q = 1 TeV, 
then there are only two additional parameters. At and tan (3. We could just as easily substitute 
tan/5 for i^(Qhigh)) which is probably a better calculational input parameter since it is not 
sensitively dependent on initial conditions, thresholds, and the loop order of the Yukawa RG 
equations, etc. Of course to precisely translate a Y(((5high) back into a tan/3, one must treat the 
above carefully. We will provide approximate tan /3 values corresponding to particular choices of 

^(Qhigh)- 
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Figure 2: Evolution of Ma/oia as a function of scale, for Yt{lQ^^ GeV) = 1, and At = 
M3, 0, — M3, — 2M3 at Q = 1 TeV for each triplet of (solid, dotted, dashed) lines from top to 
bottom. As in Fig. |l|, the solid, dotted, and dashed lines correspond to a = 1,2,3 respectively. 
The normalization differs from Fig. |l| in that the ratio Ma/ua with Yukawa couplings included 
is normalized against the ratio without Yukawa couplings (identical to the curves in Fig. |l|). 

Figure § shows the effect of a nonzero Yukawa coupling (and scalar trilinear coupling) , with 
yj(10^^ GeV) = 1, corresponding to tan/? 2.1 using one-loop evolution of the top Yukawa 
coupling. [Two-loop corrections to the top Yukawa coupling typically alter the value of tan (3 by 
about -0.05 to maintain Yt{lQ^^ GeV) = 1.] The results for Ma/a a were normalized by dividing 
by the ratio Ma/ua extracted from an equivalent model with the Yukawa couplings set to zero 
(i.e., the results shown in Fig. The effects of the Yukawa coupling terms are therefore evident 
separately from the pure gaugino terms that enter at two- loop. Since it is possible for the gaugino 
masses and gauge couplings to be established experimentally without any detailed knowledge of 
the scalar trilinear couplings, the two-loop corrections in Fig. |l| are calculable and do not pose a 
fundamental theoretical uncertainty in extrapolating to the high scale (unlike the scalar trilinear 
couplings). This provides additional motivation to normalize against the curves in Fig. |l|. 

It is clear from Fig. ^ that different choices for At{l TeV) affect the two-loop running at the 
level of nearly 10%, if a Yukawa coupling is 0(1). In the cases where > and At < — 2M3, 
the scalar trilinear coupling runs to values at the high scale that are more than a factor of 10 
larger than M3(l TeV). These input parameters may pose serious problems related to fine- 
tuning (squark masses will be driven to similarly large values), and may also be untenable if a 
charge or color breaking vacuum is encountered at the weak scale. To see how large the Yukawa 
coupling and trilinear scalar coupling need to be to get a significant correction, consider simply 
the competition between the top Yukawa term C^Y^At, and the gluino mass term B^^g^M^. 
One immediately observes that the Yukawa terms suffer from suppression in the overall constants: 
B^^/C^ = f|, 4, 1 for a = 1, 2, 3 in the MSSM. If 53 and Yt are of the same order, then At ~ 4M3 
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Figure 3: Analogous to Fig. |, with ¥^{10^^ GeV) ¥^{10^^ GeV) = 1, and the choices At 
Af) = At = M3, 0, — M3, — 2M3 for each triplet of lines from top to bottom. 



to instigate corrections of the same order as those of the pure gaugino terms. 

In the case with large tan /3, both ¥b and 1^ evolve to large values, and thus all Yukawa 
coupling terms must be included. In Fig. ^ we show the ratios Ma/ua normalized as in Fig. |2|, 
with 1^(10^^ GeV) ~ 1^(10"'^^ GeV) = 1 corresponding to tan/3 ~ 55, and the choices At = A}, = 
At = M3, 0, — M3, — 2M3 at 1 TeV. The size of the correction is somewhat larger than in Fig. |2|, 
up to nearly 15%, mainly because all three Yukawa terms are constructively contributing to the 
/3-function. 

Up to now we have been considering the evolution assuming the initial condition Eq. (|^). 
Although the simplest supergravity models and gauge-mediated models suggest such a relation 
should hold near the weak scale, it is nevertheless prudent to study other alternatives (four quite 
different examples of models with that are not expected to have "unified" gaugino masses can be 
, |2^, ^). Consider a model with gaugino masses = kaMa, such that 



found in Refs. ^1 



fciMi , k2M2 



ai 



02 



as 



(3) 



while Eq. (|) holds for Ma- In essence, we are considering two distinct models with different 
boundary conditions on the gaugino masses at the weak scale. The scaling factors ka relate a 
model that does respect the initial condition, Eq. (|l]), to one that does not. If = kaMa is 
to hold for all scales, then the scaling factors must depend on the renormalization scale. The 
renormalization group equation for the scaling factors can be obtained from the RG equations 
for the gaugino masses and gauge couplings, 



di' 



kn 
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. b 



ka) + E CI 



Ma 



Ax 



ka]tv{¥^Yx] 



(4) 
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which represents the evolution of the fraction 



Ma/aa 



Ma 



(5) 



If we take only one ka to be not equal to one, and approximate A'^/Ax ~ 1, Eq. (^) becomes 



dt 



(1 - ka] 



(167r2)2 



b; bj^a <^ x « 



tr(yjy,) 



(6) 



Technically, this is only an approximate RG equation that works well for ki or ^2, and to a lesser 
extent ^3. This is because RG equations for the scalar trilinear couplings depend on the gaugino 
masses, and we have assumed the scalar trilinear couplings are identical for both models at all 
scales to construct Eq. (|6|). However, typically the correction is not very large, and thus Eq. (^) 
is a reasonable approximation to the "true" ratio {M'^^/ Ua) / {Ma/ aa)- 

By comparing the two models using the scaling factors ka, some analytical insight can be 
gained into the effect of the nonstandard initial condition in Eq. (^) on the resulting high scale 
values. The RG equation, Eq. (|6|), indicates that ka will tend toward 1 with increasing scale if 
^ > or if simply the pure gaugino terms (positively) dominate over the Yukawa terms. To 
gain further insight into Eq. (|6|), we can approximately solve for ka by substituting "average" 
values for the quantities dependent on the renormalization scale, giving the constant which 
we schematically write as 



Then the solution to Eq. (tel) becomes 



(7) 



ka{Qf) 



1 + {kaiQi) - 1) 



(8) 



With sufficiently small Qi/Qf and Z^ > 0, it is obvious that ka converges to one as the scale is 
increased. Alternatively, with sufficiently small Qi/Qf and Z^ < 0, ka will diverge either toward 
very small or negative values (if initially /cq < 1) or very large values (if initially ka > 1). The 
mundane, but more likely possibility is that 



Qi_ 



(9) 



and so the initial conditions in Eq. (^) are preserved under RG evolution. For Qi/Q f = 10 
we find that Z^ must be 



Zk > (20,100) for 




< (0.95,0.8) , 



(10) 



10 



1.05 



1.00 




0.95 



0.90 



0.85 



io' lo' io' io' io' lo" io''io'°io'' io''io''io"io''io 
scale [GeV] 



,16 



Figure 4: The ratio Ma/ da is plotted versus scale, for a model with tan /3 = 5 and (— Mi/ai) = 
M2I0L1 = Ms/as at 1 TeV. The curves are normalized by dividing by a model with the identical 
initial conditions, except the sign of M\ is positive. The set of solid lines correspond to the 
normalized (-Mi/ai), with A^{\ TeV) = 2IM3I, IM3I, 0, -jMsl, -2|M3| from bottom to top. 
The dotted line corresponds to the normalized ratios M2/oi2 and M^/a^, which are independent 
of A^. 

and thus significant deviations to the initial conditions are not expected after running from the 
weak scale to the high scale. We have also confirmed this using numerical calculations. However, 
if Zfc were very large and positive (a numerical value of several hundred), then ka will rapidly 
evolve toward 1. To obtain a very large Zj,^ one is faced with large two-loop corrections to the 
RG equations. To be sure perturbation theory is still valid, one should go to higher orders and 
check that the three-loop terms are indeed smaller than the two- loop terms. In Ref. |2^] this was 
done for models with semi-perturbative unification, finding the interesting result that the initial 
condition on the gaugino mass ratios is not maintained under RG evolution, consistent with the 
above discussion. 

Another possibility that can be examined with the formalism above is to allow Mi and/or 
M3 to be negative. This is possible because in general there can be nonzero phases associated 
with Ml and M3 (the phase of M2 can be chosen to vanish). We limit ourselves to gaugino 
masses that are real, so that we can concentrate on merely a few different sign possibilities. To 
examine the impact of a particular sign choice on the RG evolution, assume the initial condition 



where we will examine the two cases Mi < 0, M2,M3 > and M3 < 0, Mi,M2 > 0. In Fig. § 
we plot the ratio of Ma/ua in a model with Mi < normalized to an equivalent model with 
Ml > 0, and all other initial conditions the same (we take the intermediate value tan (3 = 5 for the 




(11) 



11 




°-^°io' 10^ 10= 10' 10' lo'^ lo'io'io^Mo'io'^io^io^'io^' 
scale [GeV] 

Figure 5: Analogous to Fig. |, with Mi/ai = M2/a2 = (-Ms/as) at 1 TeV. The set of 
sohd, dotted, and dashed hues correspond to the normahzed Mi/ai, M2/a2, (—Ms/as), with 
A^{1 TeV) = 2IM3I, IM3I, 0, -IM3I, -2IM3I from top to bottom for (-^^3/03), and from bottom 
to top for the other ratios. 

purposes of this example). Since the dominant term in the evolution is g^M^, it is not surprising 
that M2I0L2 and M3/a3 are virtually unaffected by the sign choice of M\. However, M\la\ is 
significantly affected. Setting k\{\ TeV) = —1, one finds that Eq. (^) captures essentially all 
of the difference between a model with and without a negative M\. The size of the effect after 
running 13 orders of magnitude in scale varies depending on the competition between the pure 
gauge terms and the Yukawa terms in the RG equation. In the case where < 0, for example, 
there is a partial cancellation between the pure gaugino terms (dominated by 53M3 > 0) and the 
Yukawa terms, thus the RG evolution is not as pronounced. 

In the case where M3 < and Mi,M2 > 0, and the initial condition Eq. (|Tl| ) holds, there 
is a more dramatic effect on the RG evolution of all the gaugino masses due to the dominance 
of the g'^M^ term. In Fig. ^ we illustrate the RG evolution for the same initial conditions as in 
Fig. ^ but with M3 < instead of Mi < 0. We observe that the evolution of (— M3/a3) is not 
as significant as the evolution of (— Mi/ai) (for this example), mainly because the g'^M^ term is 
absent in Eq. (|6|) for a = 3. (An equivalent way of thinking about this is that the g^M^ dominates 
the RG evolution for Al-^/a^, which implies an approximate sign independence, see Eq. (|4^).) 
This also implies the Yukawa terms dominate over the pure gaugino terms in Eq. (^), giving the 
logarithmically increasing but well-defined separation between the different (— M3/a3) curves as 
Ax is increased. However, we should point out that Eq. (P) only approximately accounts for the 
RG evolution of the ratio; this is because, as noted above, Eq. (^) does not account for the Ma 
sign dependent part of the evolution of Ax- 
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2.2 Weak scale thresholds 



We have up to now avoided the issue of weak scale threshold corrections by judiciously 
choosing the initial conditions, including the gaugino masses, scalar trilinear couplings, etc., to 
hold at a scale above all weak scale thresholds. The motivation for this is that one can always 
match a weak scale theory including full threshold corrections to the initial conditions we gave 
in the examples above. In practice, this matching can become rather complicated |34, ^ 
if one demands a high level of precision. However, it is important to understand the origin of 
the uncertainties associated with weak scale thresholds, as well as recognizing that, for example, 
measured (pole) masses must be translated into renormalized masses, and the corrections can be 
large (especially for the gluino |2C, 34, |3^). Note that we have implicitly assumed the running 
gaugino masses are specified in the scheme appropriate for supersymmetry, namely dimensional 
reduction with modified minimal subtraction (DR) |22|. 

Threshold corrections arise from decoupling heavy particles from the spectrum by creating 
an effective theory without the heavy degrees of freedom that matches near the scale of the heavy 
particles. In weak scale supersymmetric theories one is interested in decoupling sparticles that 
are heavier than Mz, the customary choice of scale for starting RG evolution since the gauge 
couplings are very well measured by LEP experiments. There are two methods for handling 
threshold corrections: (1) "decoupling method", where heavy particles are decoupled by altering 
the /3-functions at the mass scale of the particles, and (2) "nondecoupling method", where fully 
supersymmetric /^-functions are retained and the effects of heavy particles are "resummed" as 
corrections in the conversion from the measured to running quantity. Strictly speaking, the 
second method is best since full logarithmic and finite corrections with arbitrary mixing can be 
incorporated (although there are sometimes ways of re-incorporating finite corrections using the 
decoupling method). However, explicitly resumming logarithms using a nondecoupling method 
is really only useful when the scale of the decoupled particles is near the scale of the conversion 
from measured to running quantity, otherwise the first method should be used. For completeness, 
we give formulae to compute the one-loop /^-function decoupling for the weak scale threshold 
corrections in Appendix A. 2 . 

It is useful to first recall that the RG equation for the ratio Ma/aa is independent of one- 
loop threshold corrections implemented as changes in the one-loop /3-function coefficients. This 
is obvious from the RG equation, Eq. (^), since it is independent of Threshold corrections 
at two-loop would, however, change the coefficients for the ratio's RG equation. But, weak scale 
threshold corrections at two-loop are negligible when we are working with two-loop RG evolution 
equations [^], since no large logarithm develops (precisely because we demand the scale of the 
superpartner masses be near the weak scale). There are residual indirect effects from one-loop 
thresholds that arise because of the two-loop terms proportional to the gauge couplings, but in 
virtually all instances such corrections can be neglected. 



Take the simplest case, the gluino mass |20, 36 1. It is easy to show that in the limit 
ruq <^ nig <C rriq the logarithmic corrections scale with logm?/m| when the gluino (running) 
mass is evaluated at a scale Q = nig. The coefficients of the log terms precisely match the 
shifts in the one-loop /3-functions (see Appendix [A.2|). However, nonzero quark masses bring 
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additional corrections, which can only be partially taken into account using a decoupling method. 
Furthermore, the corrections to the weak gaugino masses are complicated due to the mixings 
inherent in the resultant charginos and neutralinos [ 34 , 3^ . Luckily, incorporating the logarithmic 



corrections to the gluino mass, and incorporating approximate logarithmic corrections (neglecting 
mixings) in the weak gaugino masses is usually sufficient for most purposes. Since we specified 
running masses throughout the previous discussions, it was not necessary to explicitly calculate 
the corrections from translating the pole mass into the running DR mass. However, experiments 
will ultimately measure the pole mass, and so these finite corrections must be taken into account. 
In general, extracting parameters from experiment will require a careful analysis by incorporating 
both finite corrections and logarithmic corrections (probably using a nondecoupling method as 
in Ref. fs! 



2.3 Superoblique corrections 

There is further class of weak scale corrections ||3^, |3^, 43, ^] that result when 

supermultiplets are widely split in mass. These so-called "superoblique" corrections have much in 
common with oblique corrections of the standard model [^, particularly because they do not 
decouple for scales smaller than the heaviest sparticle mass. The manifestation of superoblique 
corrections is the violation the relation g = g, where g is the coupling of gauge bosons to fermions 
and scalars, and g is the coupling of gauginos to a fermion and its scalar partner. Ordinarily 
supersymmetry enforces this relation to all orders (in a dimensional reduction scheme [^]), but 
since supersymmetry must be broken near the weak scale, differing corrections to g and g are 
expected. 

Differentiating between the gauge and gaugino coupling suddenly begs the question of which 
one ought to be used in the ratio Ma/oa- To address this, there are four cases to be distinguished. 
To simplify the discussion, consider the effects of a single supermultiplet consisting of a scalar 
q and its fermion partner q with masses rriq and rUg respectively. (The extension to multiple 
multiplets and specialization to a particular gaugino is trivial.) The four cases are (1) mq,mq <^ 
M, (2) M <C rrig, rUg, (3) nig <C M <C nig, and (4) nig <^ M <^ nig, where M is the gaugino mass. 
The first case is trivial; it is easy to show that there is no distinction between g and g when all 
matter multiplets have masses well below that of the gaugino mass. The second case is nontrivial. 



and has been studied particularly for the messenger sector of gauge- mediated models 44 1. 
Even in cases when the lightest messenger scalar is 1/10 the mass of the messenger fermion, 
one finds the fractional difference {g — g)/g ^ 10^^. It is really only the case where either a 
fermion is very heavy and scalars are light, or scalars are heavy and the fermion is light, when 
superoblique corrections are significant. Motivated by the current experimental bounds on the 
masses of the scalar partners that require nig <C nig for the MSSM matter multiplet^, we will 
focus on the latter possibility. 

The one-loop corrections to the self-energy of the gaugino are shown in Fig. ^. We are 
interested in the sfermion-fermion one-loop corrections that contain the superoblique corrections. 

^With the possible exception of the scalar partner to the top quark. 
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(a) (b) 

Figure 6: One-loop self-energy corrections to the gaugino mass from (a) sfermion-fermion con- 
tributions, and (b) gauge boson contributions. 

The self-energy due to the sfermion-fermion diagram in Fig. ^(a) is 

2g'^S{q), f\ ,^ xm2 + (1 - x)m? - x(l - 

-p dx[l — x)m 5 , (12) 



IGtt^ Jo /U^ 

where p in the momentum of the gaugino, /x is the renormalization scale, and we have performed 
DR subtraction of the infinite piece. The Dynkin index S{q) is defined by S{q)S°'^ = tTq{t°'t^) in 
a normalization where S{q) = ^ for the fundamental of SU{N), and S{q) = |(y/2)^ for U{1)y 
in the GUT normalization. The correction to the running gaugino mass due to the superoblique 
corrections can be written as 

M(M) = M(M)o (^l + £^Ag^ . (13) 

where M{M) and M(M)o are the running DR gaugino masses with (without) one-loop super- 
oblique corrections applied. The function Ag is 



,"1 xm^ + (1 - x)m? - x(l - x)[M(M)ol^ 

Aq = 2Siq) / dx{l- x) In ^ ^ ^ ^ / ^ ^ ^ , (14) 

'o At'' 



approximating ^ = M^"^^ ~ M(M)o. The above expression agrees with the squark-quark one- 
loop corrections calculated in Refs. | pO| , |4^ . Note that there are additional one-loop finite 
corrections due to gauge boson loops [shown in Fig. ^(b)] that we do not present, since they 
can be absorbed into the translation between M^"^'^ and M{M) ||l^, |2^. The correction to 
the gaugino coupling arises from wave-function renormalization Z2 of the gaugino propagator 
iZ2/{'^ — M + ie), where Z^^ = 1 — —i'^{i>) after expanding the denominator of the one-loop 

df 

propagator — M — + ie), in (/i — M). The result is a correction to the gaugino coupling 

f{M) = ~gl{M)(^l + ^J2^^) ■ (15) 

using notation analogous to that in Eq. (p^). (There are additional nonlogarithmic correc- 
tions Q to Eq. (pIsD, but they decouple when M{M) <C rriq.) Thus, the ratio Malg\ is indepen- 



dent of superoblique threshold corrections. This is an important result, because it means that 
using the coupling g (that is more easily extracted from gauge boson interactions) to evolve the 
ratio Ma/ da will reintroduce one-loop corrections for scales between M < Q < rriq. Alterna- 
tively, the uncertainty in translating g into g due to superoblique corrections is a fundamental 
uncertainty in evaluating the true one-loop invariant M/g^ near the weak scale. 
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Figure 7: An example of superoblique corrections to Ma/g^, due to setting the masses of the first 
two generations of squarks and sleptons to be niq. The ratio {ga/ga)^ = {^a/gt)/{Ma/g^), where 
the denominator Malg\ is independent of superobhque corrections, but the numerator Maj g\ is 
expected to be extracted from experiment. The sohd, dotted, and dashed hues correspond to 
a = 1,2,3. We assumed the initial condition, Eq. (||), with Mi (Mi) = 100 GeV, and evaluated 
the correction at the scale Ma{Ma)- 



The size of the correction (including finite corrections) was given in Ref. |44|, and takes the 

form 

f)^ - (-) 

where the sum is over all scalars q transforming under the gauge group associated with g (or 
g) with Dynkin index S{q). Again, we are considering the case where the scalar partners are 
much heavier than the fermions. Of course we could have also considered the opposite case (see 



Refs. [41, Q), but there is motivation, for example, to set some or all of the scalar partners of 
the first two generations in the MSSM to be very heavy, of order ~ 20 TeV or so 
These models avoid flavor changing neutral current (FCNC) constraints from experiment (that 
dominantly restrict the first two generations) by pushing the mass scale of the exchanged scalars 
sufficiently high so that universality and alignment are unnecessary. In these "2-1" models the 



correction becomes [41, Q[ 

where we have assumed the scalar partners of the first two generations are degenerate with a mass 
rriq. If we further assume the initial condition, Eq. ([l|), then the size of the correction at M(M) 
can be evaluated, as shown in Fig. 0. It is clear that using the ratio Ma/g"^ (uncompensated for 
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the superoblique corrections) implies corrections of order a few to tens of percent, depending on 
the gaugino involved and the separation of scales. 

3 Two-loop effects beyond the MSSM 

Throughout the previous section various effects were discussed that could change the evo- 
lution of Mai da in the MSSM. Extensions of the MSSM could also easily give corrections to 
Ma/oia, usually through modifications to the two-loop /3-function coefficients. In the following, 
two classes of well-motivated extensions of the MSSM will be explored. The first postulates an 
extra U{1)' symmetry, using both a "minimal" U{\)' model and an model as examples. The 
second extension postulates extra matter between the weak scale and the high scale. 

3.1 Extra 

The group structure of the MSSM can be easily extended to include an extra U{1) group. 
There are many motivations for such an extension |^^, such as a solution to the /x-problem, 
string or GUT breaking to the MSSM plus an additional U{1), etc. However, enlarging the 
group structure of the MSSM does not come without its own subtleties. The most obvious 
difficulty is to assign U{\)' charge such that all the anomalies are canceled; this is, in general, 
a difficult problem. In Appendix ^ we present the conditions on the multiplicities and charges 
of matter transforming under a U{iy that cancel the anomalies. Although the U{\)y anomaly 
conditions were not explicitly stated, cancellation is guaranteed if, for example, matter is added 
in complete representations of SU{5) or is added in vectorlike pairs transforming under the 
MSSM gauge group. It was shown in Ref. that a simple model can be constructed with 
f7(l)' charges to one generation (only) of the MSSM matter that is nonanomalous, and this will 
serve as our "minimal" U{1) example. We will also examine an Eq model with three complete 
generations of 27s, which is also well known to be nonanomalous. We will not attempt to survey 
all possible extensions that include an extra U{1), but instead consider mainly just these two 
classes of models. However, the particular breaking pattern of a GUT group, such as Eq, can 
change the results depending on the scale of the breaking and whether other group structures 
exist at intermediate scales. We will briefly comment on this in Sec. |4.2| . Nevertheless, we expect 
that the results would not be significantly different if another 1/(1)' model ansatz were chosen 
that was valid up to near the high scale. 

Before we present results for particular models, it is useful to understand the origin of 
corrections to the ratio Ma/aa due to a C/(l)' extension. First, the superpotential is assumed to 
have the term 

W C YsSHiH2 (18) 

in place of 11H1H2, where the superfield 5 is a gauge singlet under SU (3)c x SU {2)lxU (l)y . This 
in itself is a modification of the MSSM that could be considered separately, but comes with its 
own set of difficulties ||5^. We will concentrate on a model with an extra U{1) that will be broken 
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when S acquires a vacuum expectation value (vev). To avoid reintroducing the /i-problem into 
these models, the term ^HiH2 is assumed to be forbidden (by a judicious assignment of C/(l)' 
charge). When S acquires a vev, an effective = Ys{S) is generated, and thus the CP-odd Higgs 
boson and all neutralinos in the MSSM will acquire mass. The particle content of this minimal 
extension includes a Z', a new gaugino X'g, a new Higgsino S, a scalar S, and a pseudoscalar. 
The pseudoscalar associated with S is eaten by the Z' , giving it mass via the Higgs mechanism. 
Three new soft terms are introduced 

-£soft C ml\S\^ + ^{M'X'sX's + h.c.) + {YsAsSHiH2 + h.c.) (19) 

(we take ms, M' , and As to be real) to give mass to the new gaugino and the uneaten scalar from 
the MSSM singlet superfield S, and a scalar trilinear coupling associated with the new Yukawa 
coupling. The neutralino mass matrix enlarges to 6 x 6, and the neutral CP-even Higgs mass 
matrix enlarges to 3 x 3. A new gauge coupling g' exists for the U{iy , giving two new couplings 
{g' and Ys) that enter the two-loop RG equations for the gauge couplings and gaugino masses. 
There is also a new one-loop invariant M' / g'"^ . For a general set of U{1)' charges Qi assigned to 
the MSSM multiplets, we have computed]^ the two-loop renormalization group equations for all 
the gauge couplings [gi, g2, gs, g') and gaugino masses (Mi, M2, M3, M') including the effects 
of the additional Yukawa coupling Ys. The results are presented in Appendix 

Unlike the MSSM gauge and Yukawa couplings, the gauge coupling g' and the Yukawa 
coupling Ys are not determined by low energy experiments, and are therefore essentially free 
parameters. (Of course there will be constraints on these parameters in particular models from 
nonobservation data.) The size of the effects of the additional 17(1)' are therefore highly model 
dependent. In practice, the two-loop corrections are not expected to be arbitrarily large unless 
g' or Ys approach a Landau pole near the high scale. However, we will only consider models in 
which g' and Ys can be treated perturbatively throughout the energy scale of interest, which will 
limit the size of the effect. 

There is an additional parameter when two or more U{1) groups are present in the low 
energy effective theory. The kinetic terms of U{1)y can mix with those of 17(1)' with a strength 
that is a priori unknown and a free parameter (again, in particular models there are constraints 

|]. Furthermore, the mixing is not RG invariant if 

triQaQb) + (20) 

the product of charges of the two groups, U{\)a x f^(l)b) summed over all chiral representations 
in the theory is nonzero. Although it would be interesting to know if there are consequences of 
kinetic mixing on the evolution of Ma/ota, we do not expect qualitatively different results, and 
therefore we do not consider this further in this paper. In particular, we will consider extra U{\) 
models that obey iriYQ) = and have their kinetic mixing set (by hand) to zero. 



''When Qa and Ma are specified for a model with an additional (7(1), the a — A elements correspond to g' and 
M'. 



from experiment) M, M, p3, 54 
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3.2 Minimal f/(l)' model 



One of the simplest extensions of the MSSM that includes an additional 17(1)' was extensively 
studied in Ref. |4^] . They proposed a superpotential with two Yukawa couplings, 

W = YuQH2u'' + YsSHiH2, (21) 

where the up-type Yukawa Yu couples only to the third generation, and the down-type Yukawas 
Yd and Y^ are absent (actually forbidden by C/(l)' charge assignments). The 17(1)' charge assign- 
ments consistent with the cancellation of anomalies are given in Appendix |B.2| . The undetermined 
parameters of the model in addition to those in the MSSM are (in our conventions): the product 
of the U{iy charge for Hi and the C/(l)' gauge coupling Qig' , the U{iy gaugino mass M', the 
Yukawa coupling I5, and the associated scalar trilinear coupling As- 

Since U{1) couplings are not asymptotically free, any additional U{1)' coupling added to the 
theory will have the greatest effect on the RG evolution of other parameters near the high scale. 
In general, it is hard to achieve any significant correction due solely to a larger g' without running 
to a Landau pole very near the high scale. In the following, we will take the initial condition in 
Eq. (Q), to hold at 1 TeV, analogous to the analyses done in Sec. ^. Furthermore, Qig' is set to 
gi{= 52) at the high scale, M' = cM^ at the weak scale where c is a constant, and for now Yukawa 
contributions to the two-loop RG equations are ignored. With these assumptions, the correction 
to MajoLa is less than 0.004c after evolving 13 orders of magnitude in scale. Even if the coupling 
Q\g' is increased to of order one, the correction is smaller than about 0.01c. Unless M' is taken 
to be considerably larger than M3 (which would be disfavored by naturalness arguments), the 
correction due to "pure ?7(1)' gaugino" contributions to two-loop running is essentially negligible. 

However, when the Yukawa couplings are nonzero, larger deviations can be obtained without 
soft parameters exceeding of order M3 near the weak scale. In Fig. ^ we show the corrections 
to MajoLa arising from a model with y5 = 1 at the high scale, and Yt finite (we set tan/3 ~ 5; 
lower or higher values of tan/3 give corrections analogous to Figs. ^ and |3|). Again, in complete 
analogy to the MSSM studied in Sec. large Yukawa couplings typically imply that the scalar 
trilinear couplings run to multi-TeV scale values near the high scale, which causes the nontrivial 
correction in Majo-a- The size of the corrections are, in order from largest to smallest, M2/a2, 
Mi/ai, M4/a4, Ma/as. The ratio Ms/as receives the smallest correction because the /3-function 
coefficient = 0, and thus it is only through Yg feeding into Yu that there is any correction at 
ah. 

3.3 //-model of E^, 

There has been an enormous amount of work studying the phenomenology of superstring 
theories, particularly the group structures that emerge from compactification. One distinct 
possibility is Eq ||5^, |5^, |5^, |6^, which can arise from under suitable breaking. E^ 
itself is rank-6, and one commonly studied breaking pattern is E^ S'O(IO) x C/(l) — > SU{5) x 
C/(l) X C/(l). Although the effects of two additional f7(l)'s could itself be a particularly interesting 
possibility to study, we restrict ourselves to considering the rank-5 77-model SU{5) x U{l)r^, and 
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Figure 8: Evolution of Ma/aa as a function of scale, for 15(10^^ GeV) = 1, and At = As = 
M3, 0, — M3, — 2M3 at Q = 1 TeV from top to bottom for each triplet of lines in (a), and for each 
line in (b). The solid, dotted, and dashed lines correspond to a = 1,2,3 respectively in (a), and 
the solid lines correspond to a = 4 in (b). The ratio Ma/aa has been normalized against an 
equivalent model with the Yukawa couplings set to zero. 



further assume SU{5) breaks near the scale where Eq was broken. (Changing the scale where 
the GUT breaks, and other possible breaking patterns will be discussed in Sec. |4.2[ ) The U{l)n 
is a mixture of the two U{l)'s in Eq that is assumed to survive]^ to near the weak scale. This 
provides a well- motivated example of a model with an extra U{1) in which the charge assignments 
and overall normalization are determined (by embedding into Eq; see Table Q for the charge 
assignments) . 

In particular, we consider an r/-model with three generations of 27s, which is the smallest 
number of representations that incorporates the MSSM and is nonanomalous. Unlike the minimal 



U{iy model considered above in Sec. 3.2, all three generations are charged under the 11(1)' , and 
Yukawa couplings for b and r are present. Although the latter is not expected to affect the results 
of evolving the ratios Ma/aa to the high scale (unless tan (5 is large), it does overcome a limitation 
of the minimal model. Three complete 27 representations of Eq include three generations of the 
following: MSSM matter (Q, u^, dP, L, e'^), down-type and up-type Higgs doublets {Hi, H2), 
two color triplets {D, D'^), a right handed neutrino {i'^), and a MSSM gauge singlet {S). The 
superpotential for the model is 

W = YtQH2U^ + YdQHi(f + YeLHie^ + YsSHiH2 + YdSDD^ . (22) 

Additional terms are possible, and indeed some are probably necessary for a phenomenologically 



viable model [32, 60 1. However, additional dimension-3 terms in the superpotential with small 



Yukawa couplings will not affect the results presented below. To prevent large tree-level flavor 



''The other f/(l) orthogonal to U{\)r) is assumed to break at the high scale. In this section, we neglect effects 
of this high scale U{V). 
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changing neutral currents, only one set of Higgs {Hi, H2, S) is assumed to couple to the quark 
and lepton superfields. To break the U{l)rj, give mass to the color triplets, and generate an 
effective ^ term, the scalar component of S is assumed to acquire a vev. While the details of the 
vacuum structure (including avoiding possible charge/color breaking minima) are very important 
to be able to construct viable models of Eq, we will not discuss this further^ Our intention is 
to take the matter content of an Eq model with the Yukawa couplings given in Eq. ([2^), and 
examine the consequences for the evolution of Ma/oa- 

The RG equations for a three generation Eq model are given in Appendix including the 
two- loop terms proportional to Yg and in the gauge coupling and gaugino mass RG equations. 
It is well known that an Eq model with three complete 27s does not preserve gauge coupling 
unification, due to the extra Higgs doublets and color triplets. Although it is possible to add 
extra matter to the theory to bring the gauge couplings back into alignment at the high scale 
(see e.g. |^, ^), we defer a discussion of the effects of extra matter (beyond that needed to 
fill a 27 of Eq) to the next section. 

There are four main effects that could alter the running of Ma/ eta in this model: the extra 
matter, the U{\)n gaugino mass, and the Yukawa couplings Ys and Yq. The simplest way to 
illustrate the effect of the extra matter is to set the gauge coupling g' and the Yukawa couplings 
Ys and Yd to be small. Although some of the extra matter will obtain masses proportional to 
these Yukawa couplings, we will simply work in the approximation where a common scale Qextra 
can be chosen for the masses of the extra matter. By this we mean the two generations of pairs of 
Higgs doublets that do not couple to quarks and leptons (the "unHiggses" ) , all three generations 
of the color triplets, the right-handed neutrinos, and the MSSM gauge singlets. In Fig. |9|(a), we 
show the evolution of Majo-a in a model where, again, the initial condition in Eq. (|l]) is assumed. 
The ratios MajoLa are normalized against an equivalent model with Qextra > -^^unif; the latter 
is equivalent to a model without the Eq matter content (at the weak scale). Since only the 
extra matter can alter the evolution (when g',Ys,Y£) are small). Fig. ^(a) shows that the three 
generations of color triplets induce the largest correction manifested in M^/a^, although sizeable 
corrections are also present for Mi/ai, M2/a2- The gauge couplings gi and 52 are nearly 30% 
larger atQ Munif ~ 3 x 10^^ GeV, while g^ is nearly 70% larger, as compared to a model without 
the extra matter. In practice this means that the effects of the additional f7(l),, coupling g' 
and the Yukawa couplings Ys, Id are not nearly as significant to the RG evolution of Ma/oa, 
since they compete against the considerably larger gauge couplings. If we include the effects 
of the U{l)r] gaugino, the shift in the curves of Fig. ^a) is at the level of about 0.007c, where 
M' = cM^ at the weak scale. Including Yukawa couplings Ys{Munii) ~ ^^(-^unif) ~ 1) with 
scalar trilinear couplings in the range — 2M3 < Ax < M3, there are more significant effects, as 
shown in Fig. ^(b). Again, the combined effects of large gauge couplings with larger Yukawa 
couplings push the scalar trilinear couplings to multi-TeV values at Munif. 

Up to now we have tentatively placed the extra matter at Qextra = 1 TeV. If the scale of the 
extra matter were increased to 2, 10 TeV, one finds the shift in Ma/ da at Q = Mu^if is roughly 

''The mapping of our results onto completely viable weak scale models derived from is beyond the scope of 
this paper. 

'^Munif was somewhat arbitrarily chosen to be the where gi and g2 meet, which is lower in this model. 
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Figure 9: Evolution of Ma/aa as a function of scale for an Eq model with (a) the Yukawa 
couplings Ys and Yd set to zero, and (b) the Yukawa couplings Is'(-^^unif) ~ YD(Afunif) ~ 1- In 
graph (a) the solid, dotted, and dashed lines correspond to a = 1,2,3. In graph (b) the four 
solid lines in the upper and lower part of the graph correspond to = M3,0, — M3, — 2M3 at 
Q = 1 TeV for Ms/as and M^ja-}, respectively. For comparison we also plotted Ms/as and 
-^2/02 from Fig. ^(a) as dashed lines. The ratios have been normalized against a model without 
the extra matter needed to fill complete 27 representations of E^. 

(-0.002,-0.002,-0.011), (-0.019,-0.019,-0.046) for a = (1,2,3) respectively, relative to the 
results in Fig. |^(a). Hence, increasing the scale of the extra matter by even a small amount (one 
order of magnitude of twelve), reduces the correction to Ma/aa by several percent. 



3.4 Extra matter 



In the r/-model of E^ described above, it was necessary to introduce extra matter to fill 
complete 27s to ensure the cancellation of extra U{\) anomalies. Indeed, the extra matter 
dominated the correction to Ma/aa- However, extra matter can, of course, be added without 
enlarging the gauge structure of the model or even disrupting gauge coupling unification, as long 
as anomaly cancellation with the MSSM group is ensured. A recent, widely touted example of 
adding extra matter without upsetting gauge-coupling unification is the addition of 5 + 5 pairs 
(or a 10 + 10 pair) used as the messenger sector of gauge-mediated models While we 

defer going into details about gauge-mediation (particularly the supersymmetry breaking masses 
associated with the messenger matter) until Sec. ^, this does serve as one interesting starting 
point for adding extra matter. Nevertheless, for the remainder of the section the extra matter is 
assumed to have only supersymmetry preserving masses. 

It is perhaps useful to review a few recent examples of the uses of extra matter for model- 
building and phenomenology (other than gauge- mediation) . Ref. showed that adding extra 
matter could "refocus" the gauge couplings to unify at a higher scale. The motivation was 
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an attempt to bridge the discrepancy between Munif and Mstrmg; and indeed they were able 
to show that only particular kinds of extra matter were able to do the job. In Ref. a 
"semi-perturbative" model was constructed whereby the gauge couplings became large, but still 
perturbatively calculable, at the high scale. They added sufficient 5+5 and 10+10 pairs to nearly 
saturate the /3-functions, and found the interesting result that the gaugino mass ratios Ma /Mb 
at the weak scale could be quite different from the canonical expectations of models without the 
extra matter. By adding matter in complete representations of SU{5), gauge-coupling unification 
(to one-loop) is not upset. However, as pointed out in Ref. |Q, extra matter can be added in 
other combinations that also maintain (one-loop) gauge coupling unification, indicating that there 
are a variety of possibilities of extra matter that could be explored. Finally, Ref. |6S] showed 



that adding extra matter in complete 5 + 5 and 10 + 10 does not change^ the value extracted 
at the weak-scale, but could increase the unification scale. In particular, by taking the heavy 
multiplets to have a common mass at the high scale but split in mass at intermediate scales (by 
RG gauge corrections), they showed that the one- loop threshold corrections largely cancel the 
two-loop corrections in the /3-functions. Although this is important to determine (separately) the 
gauge couplings or gaugino masses at the high scale, the one-loop threshold corrections cancel 
in the ratio M/a. Thus, in the examples to follow, we will generally assume the extra matter is 
introduced with a common mass at a single intermediate scale. 

Extra matter can be safely added in vectorlike pairs, allowing an arbitrarily large super- 
symmetric mass term and ensuring the cancellation of anomalies. We therefore have two sets 
of parameters; the quantity (and type) of extra matter, and the scale(s) where it is introduced. 
Using the results in Appendix extra matter can be easily incorporated by merely shifting the 
/3-function coefficients at the scale of the new matter, creating a new effective theory. Increasing 
the amount of extra matter can be compensated by increasing the scale where the matter is 
introduced, and so following Ref. [28| we introduce an effective amount of extra matter nefj that 



(in particular) need not be an integer. In this way, extra matter need only be added at one scale 
(that we take to be near the weak scale ~ 1 TeV), and then the size of the deviation in Ma/ota 
is computed as a function of the continuous parameter neff. 

The type of extra matter introduced could be of many different varieties, however there 
are constraints. If we restrict ourselves to vectorlike matter, then there could be vectorlike 
pairs of MSSM matter (e.g. Q + Q), or limitless varieties of exotics that may (or may not) fill 
representations of common GUT groups, and that may (or may not) preserve gauge coupling 
unification. We will consider only four scenarios of adding extra matter, and we decided to 
examine scenarios that do not upset gauge coupling unification (at least to one- loop). The 
scenarios are: 

(23) 
(24) 
(25) 
(26) 



Set 1 : 




= n5 + rig 


Set 2 : 


"-eff 


= nio + njQ 


Set 3 : 




= 2nQ = 2nd = 4ne 


Set4 : 


"-eflf 


= 2nQ = Auu = 2nL ■ 


*With up to intermediate values of the 


gauge 


couplings at the unification scale 


expectations are quite different [ESL [70| . 
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Figure 10: Ma/aa evaluated at Q = 10^^ GeV with four classes of extra matter, normalized to 
a model without extra matter. The solid, dotted, and dashed lines correspond to a = 1,2,3. 
Graph (a) contains ngfr number of 5 + 5 (triplet of lines on the right) and nefj number of 10 + 10 
(triplet of lines on the left), and graph (b) contains riefr number of Set 3 (triplet of lines on the 
right) and n^s number of Set 4 (triplet of lines on the left). 



Set 1 and Set 2 are the common extensions of adding vectorlike matter in complete SU (5) 
representations 5 + 5 and 10 + 10 pairs, respectively. We also consider adding 

Set 3: (3, 2, i) + (3, l,i) + 2 X (1,1,1) + conj. (27) 

and 

Set 4: (3,2,i) + 2x (3,1, -|) + (1,2, -i) + conj. (28) 

which do not apparently form complete representations of any simple GUT group (Set 3 was first 
considered in Ref. [168| ). Nevertheless, Sets 3 and 4 ensure the one- loop /3- functions are shifted by 
the same amount independent of the particular coupling, and therefore one-loop gauge coupling 
unification is preserved. 

In Fig. 1^ we show the effect of adding number of pairs of the additional matter in 
one of the above scenarios. The curves are cut off when any coupling a exceeds the somewhat 



arbitrary value 0.5; beyond this value, perturbative calculations are probably not reliable [28|. 
It is clear that the size of the correction to Ma/aa by adding extra matter can be very large, 
in comparison to all the previous effects that have been considered. However, this occurs only 
when the gauge couplings are near their semi-perturbative value, and indeed we were unable to 
find corrections much larger than 20% to 30% without a gauge coupling exceeding a ~ 0.2 at 
10^^ GeV. In all cases above, the largest correction occurred for M^/a^ due to our choice of 
adding extra matter that preserves one-loop gauge coupling unification, that implies remains 
the largest gauge coupling throughout the RG evolution. 
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The effect of increasing the mass Mextra of the extra matter can be easily approximated by 
rescahng by a factor 



Qhigh 

"eff ^ -^^n,s, (29) 

-^^cxtra 

where Qhigh ~ 10^^ GeV. If the amount of running is reduced Qhigh < 10^^ GeV, the rescahng 
of rics is more comphcated because the three gauge couphngs are no longer expected to be equal 
at the scale where the running is stopped. 



4 High scale corrections 

There are a variety of effects that can upset the ratio Ma/ eta near the high scale. In the 
following, we will examine the well-known effects of GUT scale threshold corrections [^] , discuss 
issues related to the scale of GUT breaking, examine Planck scale corrections, and finally discuss 
effects unique to supergravity. 



4.1 GUT scale thresholds 

If gauge coupling unification is indeed indicative of a grand unified theory, then many ad- 
ditional heavy fields near the high scale will be present. In SU{f>), for example, there are 
the heavy gauge boson and gaugino remnants of the adjoint, transforming as (3,2,|) under 
SU{?,)c X SU{2)l X f7(l)y, as weh as heavy Higgs multiplets (from, e.g. the 24 and 5 -|- 5 of 
Higgs), and possibly heavy chiral multiplets as well. The corrections to the ratio Ma/oia due to 
integrating out the heavy remnants of GUT breaking were calculated in Ref. ||7l[, and we will 
follow their discussion closely. There are two classes of corrections to the ratio Majcua'- wave 
function renormalization and mass renormalization. For a heavy chiral or vector supermulti- 
plet with a mass Mheavyj the wave function renormalization corrections give rise to logarithmic 
corrections to Ma and Qa that scale as 

~£ln^, (30) 
47r g^uT 

after integrating out the heavy fields near the GUT scale. Of course this is completely equiva- 
lent to modifying the /3-function coefficients at the scale of the heavy fields, and therefore the 
logarithmic corrections cancel in the ratio Ma/ota |71, p^] to one-loop. (Two-loop corrections 



can be neglected since we are assuming the heavy fields are near the GUT scale, and so no large 
logarithms are present.) 

The second class of corrections arise when soft supersymmetry breaking masses are intro- 
duced in the mass matrices for the heavy fields, leading to mass renormalization. In a general 
supersymmetric theory with soft breaking, we would expect soft terms for all sparticles with a 
scale of order the breaking scale MsusY; presumably near the weak scale. The soft breaking 
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masses are ordinarily insignificant compared with the supersymmetry preserving masses induced 
from GUT breaking. However, the corrections to Ma/aa take the form W% 



Ma ^ ^Wgut) ^ 1 / 2C(G.)[M5(Qgut) -5m] + j: S4R)Bn] (31) 

where M5(Qgut) and a5((5GUT) are the unified gaugino mass and gauge couphng at the GUT 
scale QguT) ^itT' is the mass of the fermion component of the Nambu-Goldstone multiplet in- 
duced by supersymmetry breaking, and Bji are the standard i?-terms in the soft supersymmetry 
breaking Lagrangian for the chiral representations R. The group theory factors are the quadratic 
Casimir C{Ga) for the adjoint Ga, and the Dynkin index Sa{R) for the chiral representation R, 
both for the group associated with the "a" gaugino. When the supersymmetry breaking masses 
M5, 6m and Bji are roughly the same scale as the gaugino masses, the correction is typically of 
order a few percent. 

In the case where the GUT scale chiral multiplets fill up a complete representation of the 
GUT group, the corrections J2a ^aiR)BR are finite but universal. This is completely analogous 
to integrating out the messenger sector of gauge- mediated models, and indeed the calculation 



of the gaugino mass in gauge-mediated models ||65|, [Gg] is identical to the calculation above, in 



the limit that the mass splitting of the scalar components of the messenger superfields is much 
smaller than the associated fermion mass. The presence of heavy fields with supersymmetry 
breaking masses is a bona fide possibility in itself, and the heavy fields need not be associated 
with a GUT nor do they need to have masses near the GUT scale (although a GUT does provide 
a strong motivation for the existence of such fields). Thus GUT scale corrections to Ma/aa 
are indeed possible, but the size of the effect depends solely on the size of the supersymmetry 
breaking masses associated with the heavy fields. 



4.2 Breaking GUTs at other scales 

When we discussed the Eq GUT in Sec. js!^, it was assumed that the breaking pattern 
Eq SU{5) X U{1) X C/(l) ^ SU{3)c X SU{2)l x [/(l)y x [/(l)^ occurred near the high 
scale. Extrapolating Ma/aa (for 1 < a < 4) was therefore well-defined up to this GUT breaking 
scale. However, it is very possible that the breaking of the GUT group does not go directly to 
the standard model (plus extra [/(l)'s), but instead occurs in stages. For example, the rank-6 
group Ed could break to either 50(10) x U{1) or [SU{?,)f, the latter breaking to 5C/(3)c x 
SU{2)l X SU{2)r X U{1)b-l and then to the MSSM group. The rank-5 group 50(10) could 
break to SU(5) x U{1) or 5C/(4)c x SU{2)l x SU{2)ji, and then there are various permutations 
of breaking the latter 5C/(4)c and SU{2)r. The phenomenology of these breaking patterns has 
been well-studied (for a review, see Ref. 0]), and for our purposes it is important to recognize 
that the ratios Ma/aa may only have meaning up to some intermediate scale, where part or 
all of the MSSM group is embedded in a more complicated (possibly higher rank) structure. It 



was pointed out in Ref. |73] that despite the possibly complicated breaking patterns in 50(10), 
the expectations for the gaugino mass ratios at the weak scale to one-loop were identical to 
the case where the breaking is directly to the MSSM group. Two-loop corrections are probably 
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small, unless the gauge or Yukawa couplings become large, or unless the two-loop coefficients 
are significantly altered by the extra matter in higher dimensional representations needed to 
accomplish the various stages of breaking. Above the breaking scale, the ratios Ma/ua are 
subsumed into "new" one-loop invariants M/a associated with the group at that scale [e.g. 
Ms/as for SU{5) or Mio/aio for SO{W)]. Although we have not attempted to analyze the 
various (two- loop and other) effects that could disrupt these GUT one- loop invariants, there 
would be nonzero corrections. 



4.3 Planck scale corrections 



It is well known that physics near the Planck scale can affect low energy predictions of, 
for example, the gauge couplings through higher dimension operators |7^ , [75| ]. These analyses 
presume quantum gravity induces nonrenormalizable (dimension > 4) terms suppressed by the 
Planck scale^, that can take the form |74] 



^ MWWa^) + h.C.]p , (32) 

where W"' and S are gauge and Higgs superfields and the trace is over the gauge group generators. 
If S acquires a vev for its scalar component S (through, e.g., the breaking of a GUT gauge group), 
then the kinetic terms for the gauge fields are modified 



^Ainetic = -^tr(F^,F^''^5) , (33) 



where F^jy = F^y^'^ is the field strength tensor of the unified gauge group. To bring the kinetic 
terms back into a canonical form, the fields are rescaled, and therefore the gauge couplings must 
also be rescaled. This is the well-known Planck scale correction to the gauge couplings, that has 
been discussed in Refs. ^ |7|, |, |7|, ^, |7|1. 



The Planck scale operator not only gives modifications to the gauge kinetic terms, but also 
modifies the gaugino kinetic terms (and auxiliary fields) 79, 80 1 



SC^nctic = ^tr{iXTpXS + \D^S) , (34) 

where A is the 2-component gaugino field, D is the auxiliary field, and T)^ is the covariant 
derivative. To bring these terms back into a canonical form, the gaugino (and auxiliary) fields 
must also be rescaled. Hence, the soft term ^M{\\ + h.c.) is rescaled, giving a correction to the 
gaugino mass simultaneously with the correction to the gauge coupling. It is simple to show that 
the correction to M cancels in the ratio with , and therefore Planck scale operators of the form 
in Eq. (^) do not affect the ratio M/a. This is easily generalized to Ma/ota- Furthermore, any 
rescaling of the gauge kinetic terms (from higher dimension operators, or any other source) must 
be accompanied by rescaling of the gaugino kinetic terms by the same amount and therefore 
cancels in the ratio M/a, if supersymmetry is to be preserved. 

^Normally taken to be the "reduced" Planck scale M = Mpi/\f^ 
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4.4 Supergravity effects 



In the previous section concerning Planck scale corrections of the type in Eq. (|32D, the 
operator was implicitly assumed to be globally supersymmetric. When this is generalized to a 
locally supersymmetric theory (supergravity), the Lagrangian becomes more complicated 
Specifically, the gauge kinetic terms plus gaugino soft terms take the form 

C C -iRefabF^^F'^'^' - iRe/„feA>A^ + le-^'/^G^G-y.^TA' (35) 

where two new functions are introduced; the gauge kinetic function fab, and the Kahler potential 
G, with derivatives denoted by G* = dG/d(/)i and Gj = d'^G/d(j)*'^d(i)j. (A is the 4-component 
Majorana gaugino field and (pi is the scalar component of the chiral superfields $j in the theory.) 
The gauge kinetic function fab can be thought of as parameterizing the dimension > 4 Planck 
scale operators. The Kahler potential is indigenous to supergravity, and is a (real) function of 
the scalar components of the chiral superfields in the theory. Since its form is a priori unknown, 
the coefficient of the gaugino soft mass term is unknown. Scalar fields with vevs near the Planck 
scale are expected to modify the Kahler potential, and therefore modify the gaugino masses 
independent of the kinetic terms |8^, |7^. Thus, the relation Ma/aa is expected to be modified 
when embedded in a supergravity theory with a non-minimal Kahler potential. This should not 
be surprising, since universal nonzero gaugino masses are accomplished by specifying fab = SSab 
and Gj = —dj (a flat Kahler manifold), and determined by the gravitino mass proportional to 
exp(— (G)/2). More detailed computations of the gaugino masses in stringy supergravity models 



were carried out in Ref. 1821 to which we refer the interested reader. 



5 Gauge-mediated masses 

The final class of "corrections" to Ma/aa that we will consider is from extra matter with 
supersymmetry breaking masses. Ordinarily supersymmetry breaking is assumed to occur in 
a "hidden" sector, which is assumed to contain matter with no MSSM gauge charges. Conse- 
quently, only gravitational interactions exist between the hidden sector and the MSSM, and it 
is through gravity that supersymmetry breaking is communicated to the "visible" sector; this is 
the canonical "supergravity-inspired" model. An alternative possibility is that supersymmetry is 
broken in a sector that does have couplings to the MSSM, and thus ordinary gauge interactions 
can communicate supersymmetry breaking to the MSSM. This scenario has been the study of 



an enormous amount of recent work 165, 36 



In gauge-mediated models, a "messenger" sector is introduced that has both supersymmetry 
preserving masses and supersymmetry breaking masses. The supersymmetry preserving masses 
ensure that the messenger masses are considerably larger than the weak scale, while the super- 
symmetry breaking masses can successfully induce weak scale masses in the MSSM sparticles. 



The canonical approach |65] is to assume the messengers form complete representations of SU{5) 



to ensure gauge coupling unification is not disrupted at one-loop (however, see Ref. fB^])- The 
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Figure 11: The one- loop gaugino mass induced by heavy messengers. The crosses represent the 
relevant mass insertions. 



superpotential for the messenger sector takes the form 



W 



(36) 



where $ transform as a vector like representation of the MSSM gauge group, such as 

(3,l,-i) + (l,2,i) + conj. 



(37) 



that make up the 5 + 5 of SU{5). Supersymmetry breaking occurs when the auxiliary component 
of the gauge singlet S acquires a vev Fs, but we also assume the scalar component of S also 
acquires a vev (S) to give $ and $ (large) supersymmetry preserving masses. The spin-1/2 
components of ^ and $ are combined as one single Dirac fermion, with a mass mrif = \{S). The 
two complex scalars from <I> and $ acquire masses mi and m2, with m\ 2 = ^tt-j + AF^. In general, 

the effect of such terms is 

2 I 0™2 A^2 



additional (higher dimensional) interactions could be added to Eq. (| 

to modify the relation between mi^2 and m/, parameterized by Str M^^gg = 2mf +2m2— 4mj [33| 

The one-loop messenger contributions to the gaugino masses have been presented in Refs. p5| , 
|6^ ], and generalized to the case Str M^g^g 7^ in Ref. The relevant diagram is shown in 



Fig. 11, which is equivalent to the mass renormalization diagram due to heavy chiral multiplets 



discussed in Sec. and is the supersymmetry breaking analog to Fig. [^(a). The result for the 
self-energy can be written as 



a 

47r' 



-m 



/5($)sin20 X 



y2 



1-^2 



In 2/2 



yi 



1-2/1 



In yi 



(38) 



where yi = mf/{p^ + m^), mj is the messenger fermion mass, mi and 771-2 are the messenger 
scalar masses, 6 is the angle that diagonalizes the mass matrixf^, and S'($) is the usual Dynkin 
index for the messenger superfield The correction to the gaugino mass can be approximated 
as simply 



M = Mo + S(Mo 



(39) 



where M and Mq is the gaugino mass with and without the supersymmetry breaking sectors 
integrated out. The expression given in Eqs. (^), (^) is a trivial generalization of the formulae 



given in Refs. [68, S3|, and reduces to the latter formulae in the limit that the gaugino mass is 
induced entirely from messenger interactions. In fact, the p"^ correction in the denominator of yi 
is negligible when the messenger fermion mass is much larger than the gaugino mass, which is 
the usual situation. 



We assume the scalar mass matrix is real. 
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The size of the correction to the gaugino mass is dependent entirely on the size of the 
supersymmetry breaking mass \/Ts in the messenger sector. In the case of heavy chiral multiplets 
with supersymmetry breaking masses of order the weak scale [^], the correction is of order a 
few percent. This is because the one-loop self-energy correction S(Mo) is suppressed by a/(47r). 
In gauge-mediated models, this suppression is overcome by specifying a supersymmetry breaking 
mass of order the weak scale times iir/a, which is of order tens of TeV. Note that two-loop 



corrections to the gaugino masses have also been recently computed in Ref. |84|, and can give 
corrections up to 10% to the one-loop result if the messenger scale is of order the supersymmetry 
breaking mass. 

It is also worth considering the possibility that extra heavy vectorlike matter exists with 
supersymmetry breaking masse g0 that do not come in complete representations of SU{5), nor 
even preserve gauge-coupling unification. This class of matter would communicate its supersym- 
metry breaking to only those sparticles charged under the appropriate groups to which the extra 
matter transforms nontrivially. There are, however, some important comments. The size of the 
contribution to any given gaugino mass (and thus the ratio Ma/aa) is determined by the size 
of the supersymmetry breaking mass (and the multiplicity of the extra matter). The correction 
to the gauge couplings (at the high scale) due to the extra matter is determined solely by the 
multiplicity and the scale of the extra matter, i.e. the supersymmetry preserving massesf^. Thus, 
some classes (or multiplicities) of extra matter could be excluded if the gauge couplings are re- 
quired to be perturbative to the high scale. The size of the correction to the gaugino mass is, 
however, undetermined. Hence, the correction to Ma/oa could be arbitrarily large. 

There are, of course, phenomenological consequences if supersymmetry is broken at scales 
smaller than the canonical scale ~ 10^^ GeV assumed in supergravity models. In particular, a 
light gravitino is expected. The presence or absence of a light gravitino in collider experiments 
would help determine (or at least put a lower bound on) the fundamental scale of supersymmetry 
breaking (see e.g. ||8^, ^). However, it probably does not determine the supersymmetry breaking 
masses in the extra (messenger) matter because generically one expects complicated feed down 
mechanisms in the dynamical breaking sectors that ultimately separate the "true" scale \/F from 



6 Conclusions 



In this paper, we have studied the effects that can disrupt the one-loop renormalization 
group invariant, Ma/ua- This includes two-loop corrections in the MSSM, weak scale threshold 
corrections, and superoblique corrections. Two-loop corrections from extensions of the MSSM 
within a minimal U{1) model, an Eq model, and extra matter were also considered. In addition, 
high scale corrections from GUT thresholds, Planck scale operators, and supergravity effects 

^'^To avoid the introduction of additional complications, the supersymmetry breaking can be assumed to be 
introduced explicitly. Those who find this unpalatable can imagine gauge singlet interactions (plus possibly higher 
dimensional terms) analogous to Eq. (p6[). 

^■^Highly split supermultiplets will have additional corrections dependent on the size of the supersymmetry 
breaking masses. 
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were discussed. Finally, the effects from a heavy sector with super symmetry breaking masses 
that communicates via gauge-mediation to the MSSM were discussed. 

The two-loop corrections (both within the MSSM, and beyond) were, not surprisingly, found 
to be large only when a coupling becomes large near the high scale. Otherwise, the corrections 
are of order a few to tens of percent after evolving 13 orders of magnitude in scale. The "pure 
gaugino" terms were usually small except for the (/3M3 term, and in particular the extra 11(1)' 
gaugino in extended models typically gives a very small correction. The Yukawa terms can 
induce a larger correction (up to ten to perhaps tens of percent at the high scale), particularly 
if tan/3 is near one of its perturbative extremes. Typically, to obtain a larger correction the 
scalar trilinear couplings must run to several times their weak scale values at the high scale. The 
Yukawa terms, and hence the correction to Ma/ eta, can take either sign depending on the sign 
of the scalar trilinear coupling (although note that the sign of the scalar trilinear coupling is not 
RG invariant). The corrections from extra matter can also be substantial, becoming very large 
(factor of 1.5 to 2) if the gauge couplings approach their semi-perturbative values. 

We also showed that superoblique corrections cancel in the ratio Ma/g^, where ga is the 
coupling of the gaugino to the scalar and fermionic piece of a chiral superfield (that transforms 
nontrivially under the associated gauge group). However, since experiments measure the gauge 
couplings extracted from gauge boson interactions (at least for the foreseeable future), there can 
be a significant correction in translating Ma/aa into the true one-loop invariant Ma/g^, if some 
super multiplets are widely split in mass. 

Weak scale threshold corrections proportional to a (possibly large) logarithm cancel to one- 
loop in the ratio Ma/aa- There are, however, residual finite effects when translating measured 
pole masses into DR running masses, and effects from gaugino mixing. Thus, complete weak 
scale threshold corrections would be necessary for high precision. To compute the size of other 
effects considered in this paper, we specified the running masses at a common scale (1 TeV) 
above weak scale thresholds. This was mainly for computational convenience since the complete 
weak scale threshold corrections can become rather complicated; of course one can always do the 
matching of our initial conditions to a particular weak scale model if desired. 

High scale effects were also discussed. We reviewed GUT threshold corrections that have 
been calculated before, noting that the correction depends only on the soft breaking masses of the 
heavy GUT fields. GUT breaking that occurs in stages (or at a smaller scale) was also mentioned, 
implying the ratios Ma/oia may only be applicable up to a scale that could be considerably smaller 
than the unification scale. Planck scale corrections from one commonly considered operator were 
found to cancel in the ratio Ma/aa, because supersymmetry enforces equal rescalings of the 
gaugino mass and gauge coupling squared. Finally, corrections unique to supergravity were 
discussed, where they were shown to depend on the form of the gauge kinetic function and the 
Kahler potential. 

Generating gaugino masses by gauge-mediation was discussed in the final section. We calcu- 
lated the one- loop gaugino masses induced by a heavy chiral multiplet, in agreement with other 
well-known calculations when the mass is induced entirely from gauge-mediation. Since the ex- 
istence of a messenger sector and the size of the supersymmetry breaking masses is unknown, 
there can be arbitrarily large "corrections" to the ratio Ma/aa- If gauge-mediated supersym- 
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metry breaking at low energy is the only source of the gaugino mass, then the two-loop effects 
considered in this paper are expected to be much smaller than a supergravity model, because 
the RG evolution is considerably reduced. 

We should point out that many of the effects that we have discussed could occur simul- 
taneously (either constructively or destructively). Determining the origin of the corrections to 
Ma/aa would therefore be quite difficult. However, with sufficiently precise measurements of the 
weak scale spectrum it is expected that some of the possible disruptive effects could be ruled out. 
Extracting the values of Ma/ota unambiguously at the high scale relies on finding techniques that 
can set limits or exclude the myriad effects discussed in this paper. Since these effects include 
corrections from e.g. extensions of the MSSM and high scale effects, an alternative point of view 
is that the ratios Ma/aa are potentially sensitive to such physics. 

Since the supersymmetric parameters M, y, A) are sufficient to compute the RG evolu- 
tion of Ma/ eta, there was no need to specify any other parameters such as squark or slepton 
soft masses. This allowed us to compute various corrections without any significant model de- 
pendence, although particular models were used to illustrate the effects of an extra U{1). Since 
no attempt was made to construct complete extended models, the corrections to the MSSM ra- 
tios Ma/aa should be considered as generic expectations. Some particular weak scale scenarios 
(or particular sets of parameters) that we have presented may be untenable when embedded in 
a more complete model for reasons such as vacuum stability, electroweak symmetry breaking, 
constraints from colliders, and constraints from cosmology. Nevertheless, unambiguous extrapo- 
lation to the high scale can be expected to be difficult regardless of the model framework, given 
the varied effects that impact the one-loop invariant Ma/aa- Careful consideration of the effects 
discussed in this paper would be necessary to make an unambiguous connection between weak 
scale and high scale physics. 
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Appendix A: Relevant /3-functions 
A.l Gauge and gaugino /5-functions 

The /3-functions for the gauge couplings, gaugino masses, and the ratio Ma/aa in the DR 



scheme [22| have been given before in Refs. [f9, 2C, 21, 26 1 and we generally follow the notation 



32 



of Ref. !§: 
d 



dt 



9a 



167r2 ^ (167r2)2 



167r2 <^ " ' (167r2^2 
Consequently, the two-loop RG equation for the ratio Majota is 

,(2) 



(40) 
(41) 



dt gl 



(167r2)2 



Y: B)^>glM, + 5: Q-^. tr(yiy.) 

h X 



(42) 



where 1 < a, 6 < 3 for the MSSM [1 < a, 6 < 4 for the MSSM plus an additional U{1); the 
a = 4 elements of ga and Ma correspond to the gauge coupling g' and gaugino mass M' as- 
sociated with the C/(l)'], and x = u,d,e{, S, D) for the Yukawa couplings Y^, Y(i, Ye(, 1^5, Yd) 
and scalar trilinear couplings A„, A^, Ae(, A5, Ad). (The notation ^d^,. tr(yjyi,0 is understood to 
mean A^iYj^^Y^^ for the flavor indices At least the first three Yukawa couplings and tri- 

linear couplings (m, d, e) are matrices in family space, however we only consider third generation 
couplings to be nonzero. The one-loop and two-loop coefficients are written as 
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With additional matter may also come additional Yukawa couplings in the superpotential, and 
thus is potentially an 3 x m matrix for the MSSM [4 x m matrix for the MSSM plus an 
additional U{1)] with m distinct Yukawa couplings. The coefficients 5^, bab and c^x 
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in terms of the MSSM multiplets {Q,u'',d'',L,e'',Hi,H2), an MSSM gauge singlet (S), and the 
(exotic) adjoints 3 and 8. The charge assignments to SU{3)c, SU(2)l, U{1)y, and 11(1)' and 
multiphcities rii of the matter supermultiplets are shown in Table ^. To incorporate the extra 
matter that exists in a 27 of Eq, let UdQ^ nc^Q'^+noQ^+nocQ^^ and nsQ^ — > nsQ^+nycQ^^ 
in the above equations for 7^ = 0, 2, 4. Note that for the MSSM without a singlet, we could have 
ignored the distinction between ni and 77-2, and written the above in terms of just nn = ni = 712- 
However, the presence of the term YsSHiH2 as a substitute for ^HiH2 in the superpotential 
implies Qi 7^ —Q2-, and thus rii and 712 must be treated separately. 

We illustrate the utility of the above expressions by calculating the /3-function coefficients 
B^^^ and i?*-^-* for a few particular cases. First, for the MSSM 



riQ = Hu 

77-1 = n2 

713 = "8 

we obtain the very well-known result 
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Table 2: Charge assignments for the matter considered in this paper. We did not list the 3 
and 8 with multiplicities ns and ng that are in the adjoints of SU{2)l and SU{3)c respectively, 
since they are singlets under all other groups. The 77-model of Eq is a special case of a general 
U (1)' extension with the charges shown. Note that Qem = T3 + Y/2, and U {1)y is in the GUT 
normalization. 

For a model with additional 5's or 5's of SU(5) that transform as (SjQs) and (5, — Q5) under 
[SU{5), 11(1)'], one can compute the shift in the one-loop and two-loop /3-functions 
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and similarly for additional lO's or lO's of SU(5) that transform as (10,(5io) and (10, 
under [5C/(5), C/(l)'], 
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Note that we chose the simplest assignment of the 17(1)' charge, such that it is assigned vectorially 
and commutes with SU{5). Other assignments can be handled by simply taking the appropriate 
components of the decomposed SU(5) field individually, and assigning U(l)' charge as desired 
(for examples of other C/(l)' charge assignments, see Ref. [^). The [SU{3)c, SU{2)l, U{1)y] 
pieces of the above /3-functions were also calculated in e.g. Refs. ||8^, and we agree with 

their results. 

Finally, we use the above results to compute the /3-function coefficients for the r/-model of 
Eq with three 27 matter representations. The charge assignments of the U{l)rj are shown in 
Table 0. We obtain 
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(49) 
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The one-loop and two-loop /3-function coefficients were also given in Ref. [| 
completely with the results in Ref. |63[, and we agree with the results in Ref. i5 



63]. We agree 



(2) (2) (2) 

overall sign and the values of the two-loop elements i?23 > -^14 ' -^41 



except for the 
in our notation). 



A. 2 Thresholds 



In running the RG equations, we will encounter two kinds of thresholds. The first, what we 
call "fully supersymmetric thresholds" , occur when extra matter is decoupled from the spectrum 
in a fully supersymmetry way. The RG equations as given above are sufficient to remove matter 
(both the scalar and fermionic components of a chiral superfield) at any scale, provided the masses 
of the scalar and fermionic components are the same. Both one-loop and two-loop decoupling 



can be achieved; the importance of the latter for the Ma/aa ratio is discussed in Sec. |3.4 The 
second kind of thresholds are those near the weak scale, when the supersymmetric components 
of the SM spectrum are decoupled. One way to treat weak scale thresholds is to decouple 
the supersymmetric components from the /9-functions, creating non-supersymmetric effective 
theories. Expressions for the corrections to the one- loop /3-functions are given in Ref. and we 
have translated their results into the formalism we have used above: 



III —> .^ib^ -r g/tg 6 f ^ 6 eL 

^^2 ^ hf^ + l'^^l 

where / is the fermion and s is the complex scalar. In addition, the /3-function coefficients for 
the non-Abelian gauge couplings have the further changes 

62 : -6 ^ -f + 
63: -9 ^ -11 + 2% 

where - = or 1 if the renormalization scale Q is less than or greater than |M,?, - I respectively. 
Finally, note that the two-loop corrections due to weak scale thresholds could in principle be 
calculated, but we ignore these corrections since they are not enhanced by a large logarithm 



(unlike fully supersymmetric thresholds), see Ref. [37| for further discussion 
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A. 3 Yukawa and scalar trilinear /3-functions 



The general form of the one-loop RG equations for the Yukawa couplings and soft scalar 
trilinear couplings is 



dt ' 



dt 



Yj 
167r2 

1 

8^ 



E 



c Y^ 



k 



J2 CijYjAj + ^ dikolMk 
i k 



where the coefficients are 



Hk 



6 
1 

3 


13 
15 

T_ 
15 

9 



V 



15 



1 


\ 


1 1 





4 1 





1 4 


3 


2 


5 ; 


16 
3 


2(Q 


16 
3 


2(Q 





2(Q 





2(Q 


16 
3 


2(Ql 


We 


are r 



s + Qi + Q2 

D 



+ Qh) J 



(51) 



(52) 



Two-loop corrections to the 
have been 



26 



for i,j = (u,d,e, S, D) and k = 1,2,3,4. 

Yukawa couplings and associated mixings between generations. 
Yukawa couplings 0, §, 25, 2£] and the associated scalar trilinear couplings 
calculated, and expressions for the two-loop corrections to Ys and Yd could be easily calculated 
from Refs. p5| , p^ . We consistently neglect two- loop corrections to Yukawa couplings since they 
are of the same order as three-loop corrections (which we also neglect) to the main observables 
Qa, Ma, and Ma/oa relevant to this paper. 



Appendix B: Cancellation of U{iy anomalies 

B.l General conditions 

There are six anomalies associated with enlarging the MSSM gauge group to include an 
extra f7(l)' that need to be canceled if we expect a sensible low energy effective theory. These 
are the mixed gauge anomalies, [SU{3)c]^U{iy , [5C/(2)l]2c/(1)', [U{l)Y?U{iy, [U{l)fU{l)Y, 
the triangle anomaly [C/(l)']^, and the mixed gravitational anomaly C/(l)' [gravity]^. With the 
addition of the mixed gravitational anomaly, this appendix follows closely the discussion of 11(1)' 
anomalies in Ref. [^ ]. 

With the particle content, multiplicities, and U (1)' charges stated in Table |2|, the cancellation 
of the [5C/(3)e]2C/(l)', [SUi2)L?Uiiy, [C/(l)y]2C/(l)', [C/(l)f C/(l)y , [t/(l)f , C/(l)' [gravity] 2 
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anomalies require, respectively, 

= 2nQQQ + UuQu + naQd + noQn + ud'^Qd'^ 
= SuqQq + ulQl + niQi + 712(52 

= 2nQQQ + IGuuQu + ^ndQd + ^niQL + '^'^neQe + QniQi + 6n2Q2 + '^noQi) + inD'^Qo'^ 
= ngQg - 2nuQl + "-rfQd " "-lQl + rieQl - "-iQi + ''t-'^QI - ^dQe, + ud'^Q'dc 
= 6nQQ| + 3n„Q^ _^ 3„^g3 _^ 271^^1 + UeQl + 2niQ? + 2n2Q2 + ^tidQ'd + ^nocQ^c 
+ nsQ's + n^-Qi^c 

= GhqQq + SriuQu + ^ridQd + ^ulQl + UeQe + 2niQi + 2n2Q2 + SudQd + ^no'^Qo'' 
+ + n^'^Qv'^ ■ 

For completeness, the further constraint of preventing radiatively generated kinetic mixing tr(yQ) = 
is 

= uqQq - 2nuQu + ndQd - tilQi + neQe - niQi + n2Q2 - noQu + tld'^Qd'^ ■ 



B.2 Minimal model 



The model of Ref. has the superpotential 

W = YtQH2U^ + YsSHiH2 , (53) 
(devoid of Yukawa couplings to the 6 or r sector), with the MSSM spectrum and one singlet S, 



and the MSSM group plus an additional U{1)' . The authors of Ref. [49| wrote down one simple 
example of anomaly cancellation in their model, 

Qq = -^Qi 

Qu = ^{Qi — 3Q2) 

Qd = i(Qi+3Q2) 

Ql = -Q2 

Qe = Q2 — Ql 

Qs = -(Q1 + Q2), 

where the charges given above for {Q, u, d, L, e) are assigned only to the third generation (first 
and second generations have zero 11(1)' charge). To prevent radiatively generated kinetic mixing, 
we also impose the constraint 

Q2 = iQi- (54) 

Thus all of the charges are determined up to one "overall" charge, which we take to be Qi. Since 
the gauge coupling g' is always accompanied one power of some 17(1)' charge, we can therefore 
choose to rescale either Qi or g' , without loss of generality. If the 17(1)' were embedded in a 
higher rank group at the high scale, then the rescaling factor would be fixed by the particular 
embedding. 
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Note that experimental bounds on Z-Z' mixing also impose constraints on the size of g'Qi, 
and may not even allow Eq. ( [5^ ) to be satisfied in the particular scenarios where the vev of S is 
comparable to the vevs vi and V2 (from Hi and H2) [^]. The mixing angle in the Z-Z' system 
is given by 



tan2az-Z' = mI^, " ^^'^ 

We find the intriguing result that simultaneously satisfying Eq. (^4|), and zero mixing occurs 
when tan/5 = |. Of course the mixing could be nonzero but simply suppressed by a large Z' 
mass, although generally Q1Q2 > is needed for a partial cancellation in the numerator of the 
above expression [which is compatible with Eq. (^)]. In any case, for our purposes we will take 
Eq. (|5^) without specifying tan and assume that appropriate (larger) vevs of S are induced as 
necessary to avoid the experimental bounds. 



B.3 Ea model 



Using the charge assignments given in Table it is easy to show that the r/- model of Eq with 
three generations of 27s cancels all of the anomalies above. This is not at all surprising, since it 
is well known that theories with matter in complete representations of a non-Abelian group are 
nonanomalous. Three generations of 27s are sufficient to incorporate the matter content of the 
MSSM, but also include two extra generations of up-type and down-type Higgs doublets, three 
generations of right-handed neutrinos, and three generations of MSSM gauge singlets. 
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